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Single particle states in the atomic trap employing the rotating magnetic field are found using 
the full time-dependent instantaneous trapping potential. These states are compared with those of 
the effective time-averaged potential. We show that the trapping is possible when the frequency of 
the rotations exceeds some threshold. Slightly above this threshold the weakly interacting gas of 
the trapped atoms acquires the properties of a quasi- ID system in the frame rotating together with 
the field. The role of the atom-atom interaction in changing the ideal gas solution is discussed. We 
show that in the limit of large numbers of particles the rotating field can be utilized as a driving 
. . . force principally for the center of mass motion as well as for the angular momentum L = 2 normal 

' modes of the Bose condensate. A mechanism of quantum evaporation forced by the rotating field is 

analyzed. 
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I. INTRODUCTION 



The novel methods for storing atoms at very low densities and temperatures open up new opportunities 

for studying the role of the atom-atom interaction in macroscopic quantum phenomena |j] . The problems of Bose- 
Einstein condensate formation |^-|^ , the dynamical response of the condensate in the trap ^ , and the interaction of 
. the condensate with light can now be investigated experimentally. 

■ The zero dimensional geometry and the small size of the atomic traps restrict the direct observation of the most 
J — spectacular effects known from the history of superfluidity of Hell (see in, e.g., Therefore devising new practical 

methods for probing the condensate in the atomic traps becomes of crucial importance. In this regard the recent 
, suggestion [ pd| to analyze the rotational properties of the trapped atomic cloud appears to be very promising. As 
was pointed out in Refs. an analysis of quantum evaporation from the condensate can yield valuable information 
about the interatomic interaction. Therefore, adopting this analysis to the trapped gases is highly desirable. The 
very recent theoretical [|l3| and experimental analyses of a Bose-Einstein condensate undergoing variations of the 
"i ' trapping potential address the long standing question regarding the coherent versus dissipative behavior of a many 
' body system. 

O The trap jl^ where the Bose-Einstein condensation of Rb atoms was first achieved , utilizes a rapidly rotating 

magnetic field (RMF). This field, if averaged over the rotational period, creates an effective static oscillator potential 
L* ' (which is called the time orbiting potential (TOP) [p^ ) and reduces the escape of atoms from the trap due to spin-flip 
• I— I 

X 



effects |15|. Other traps Mpl (see also |16 ) do not rely on the RMF. A main assumption made about the RMF 



is that as long as the frequency of rotations ui is much larger than the frequency uio of oscillations in the TOP, the 



?H [ trapped atoms are not disturbed by the time variations of the instantaneous potential. Accordingly, the results 18 
1 obtained for the Bose-Einstein condensation in a static parabolic potential can be applied to this case as well. 

Generally speaking, the RMF should transfer energy and angular momentum to the condensate. Therefore, the 
RMF can be viewed as a possible tool for studying the dynamical response of the condensate. In this sense, addressing 
the problem of the exact description of the quantum atomic states in the trap employing the RMF, rather than relying 
on the time averaging procedure IT^ ], appears to be quite important. 

In this paper we study various aspects of the RMF: we find the exact single particle states in the trap jl^] without 
relying on the time-averaging procedure; it is shown that the RMF, if properly modulated, should excite selectively 
some modes of the condensate; we derive the Ginzburg-Gross-Pitaevskii (GGP) equation taking into account the 
effects of the quantum evaporation induced by the RMF. As an application, the decay rate of the condensate due to 
the RMF is calculated in the case of steady rapid rotations of the RMF. 

The outline of the paper is as follows. In Section II we find the exact eigenenergies and eigenstates for a single 
particle in the trap ||l5[. This solution is obtained in the frame rotating together with the RMF. The properties of 
the solution as a function of lu are analyzed. In Sec. Ill we consider the limit of large numbers of particles in the 
condensate and analyze the nondissipative interaction between the RMF and the normal modes of the condensate. In 
Sec. IV the GGP equation with dissipation due to the RMF is derived under certain approximations. The quasistatic 
solution for the rate of the centrifugal evaporation of the condensate is derived in the limit of high to. 
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II. SINGLE PARTICLE STATES IN THE ROTATING FRAME 



In our analysis of the behaviour of a single atom trapped by the magnetic field B we follow the approximation that 
the atomic spin orientation is parallel to B . Then the effective potential energy of the atom seeking the low field 
is essentially the Zeeman energy U — l/isBj, where /is stands for the Bohr magneton (we ignore the nuclear magnetic 
moment). The magnetic field of the trap [T^ ] consists of the static quadrupolar part Bg having axial symmetry 
with respect to the z-axis, and the RMF B^ rotating in the x, y plane. Representing these explicitly, one finds the 
components 

Bqx — B'yX, Bqy — B'^y, Bqz = —'2B'gZ, 

(1) 

Bbx = -BbCos{ujt), Bby = -BbSm{ujt), Bbz = 0, 

where B'^, Bb stand for the constant gradient of the quadrupolar field, and the constant amplitude of the RMF, 
respectively. Given Eq.(l), the potential energy 

C/= |MB(B, + Bb)| (2) 

depends on time. As suggested in Ref. jl^], for high uj the time dependence can be effectively averaged over, which 
results in the TOP-potential In general, the time dependence in U should result also in the nonadiabatic exchange 
of energy between the external field and the atoms in the trap. However, as will become clear from the following, for 
a time independent lu and noninteracting atoms, no such an exchange occurs between the RMF and the atoms in the 
trap. In fact, in the frame rotating together with the RMF, the single-particle Hamiltonian becomes time independent 
insuring that the atom once prepared in the pure state (in the rotating frame) will live forever in such a state. Going 
to the rotating frame implies the coordinate transformation 

x" = cos{e{t))x + sm{0{t))y 

(3) 

y' ^ - sm{9{t))x + cos{9{t))y 

where 6{t) stands for the angle between B^ and the x-direction. In the case of steady rotations, 9{t) — uut. The 
transformation (3) results in (2) rewritten in the time independent form as 



U=\fiBB'^Wix"-Xor+y'^+Az^, xo^^. (4) 
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Since the effective size of the atomic cloud is much less then Xo |15| , one can expand (4) in terms of 1/xo- This 
gives for the first two terms (linear and quadratic) 



where we have omitted the unimportant constant \xo\', the notations oj'^y — l/|a;o|, uj"^^ = 4/|a;o| are introduced, and 
the units of energy and length are employed as 

■ (6) 



° Mil' ° (M|Msi?^|)V3' 

respectively. In (6), M stands for the atomic mass. Note that in the rotating frame the stiffness of the potential along 
the ^''-coordinate is zero. This would imply that no states localized around the origin exist. However, as we will see 
below, the finite kinetic energy of the particle changes this conclusion for sufficiently large uj. 

In the rotating frame the kinetic energy acquires the Coriolis term —ujL, where L is the z component of the angular 
momentum operator and u = 9. Consequently, taking into account (5), one finds the single particle Schrodinger 
equation (?i = 1) in the rotating frame 

(7) 

= -iA + iiu{x"dy' - y'dx") + U 
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We consider first the case lu = const. Note that Eqs.(5), (7) represent a quadratic form which can be diagonahzed 
explicitly (see Appendix A). Prior to solving it let us eliminate the linear term —x" from (5). This can be accomplished 
by the transformation 

y' 1 

■0 exp(— z — )'ip(x' ,y' , z), x' = x" (8) 

which results in Eq.(7) being rewritten as 



2 „/2 



(9) 



K = --^^i^{x'dy,-y'd,,)^ 



For the case |a;| < Woy, no discrete states localized near the origin x' = y' = z' = exist in the trap. Accordingly, 
we will not analyze this case any more. For \lo\ > ujoy such states do exist. Their eigenenergies are (see Appendix A) 

Eranl ^ ^^+m- UJ^n + LOozli = i^\J 1 + '7^/2 ± 77v^2 + 77271, rj = UJoy/^^, (10) 

where m, n, I arc integer nonnegative quantum numbers, and the energy of the state with m = n = I ^ Q i& set equal 
to zero. The normalized eigenf unctions are (see Appendix A) 



Vmnlix ,y ,Z] — P 



(11) 



tl=t2=t3=0 



where we have introduced 



1 7, 1 7. *7oa; y , 



, -2 

\i{ - ii) -t'i- 2vtit2 + 2V^t3Z- 



(12) 
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+\/2[a^^ ^JuJlx' — iva^Juj^y']ti + \f2\va^JZJ{x' + ia^^ ^fuj^y'\t2- 
In Eq.(12) the parameters are 



7o = w ■ , cj2 = Vl-TWi, r] = , i^ = sign(w), 

0^2 + ti'l W 



(1-4-1^8+^)^^' 

(1 + ^ - i V8T^)l/4(l _ ^2)1/8 ■ 



(13) 



Note that the parameter a as a function of w (or rj) has the property a{—ri)a{ri) — 1. In the limit r] —^ (co oo) one 
obtains a = 1, and the eigenfunctions (11), (12) as well as the eigenenergies (10) become exactly those characterizing 
the TOP H as seen from the rotating frame. In particular, the spectrum acquires the form 

^Lni = ^oim + n) - uj{n - m) + Uozl (14) 

where lOq — '^oyl^f^ is the frequency of the oscillations in the a;, y-plane of the TOP [T^ . 

For large w, Eq. (14) is an approximation of the exact expression (10). The corrections due to the finitcncss of 
LOojbj turn out to be of the order of (wq/w)^, so that one can effectively ignore these even if uj is only a few times 
larger than lOq- Consequently, for such uj we will employ Eq.(14) instead of the exact form (10). 

In the limit \r]\ 1 from below, the solution (10)-(13) acquires features characteristic of a quasi-lD system. Indeed, 
taking this limit in Eqs. (10)-(13), one finds 



UJ+ = Ui = VScOoy, UJ- = ^UJ2^ ]J'^{l~\'n\)uJoy. (15) 
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This expression together with Eqs. (11), (12) imply that the typical extension in the y'-direction diverges as 
(1 — |f?|)~^^'* — > oo. Accordingly, the excitation spectrum (10) becomes characterized by the soft mode whose energy 
goes to zero. This implies that in this region of lo the low energy dynamical response of the system of atoms 
should exhibit the ID behavior. In this paper we will not focus on the properties of such a ID atomic gas. 

Note that the spectrum (10) and its limiting form (14) have no lower bound. For the model Hamiltonian whose 
potential is the axially symmetric TOP this is a pure consequence of the coordinate transformation (3) because 
this Hamiltonian conserves angular momentum. Accordingly, no instability with respect to a spontaneous growth 
of angular momentum L = h{n — m) (in physical units) in (14) is expected to occur. In contrast, the Hamiltonian 
(9) does not conserve angular momentum. This implies that under certain conditions such an instability could be 
realized. 

Another interesting feature of the solution (11), (12) is the common phase factor exp{i"fox'y'). Its magnitude is 
controlled by the asymmetry of the eigenfunctions in the x, y plane (in addition to the squeezing in the z-direction) . 
Also, this factor implies a very specific pattern for the velocity v = VIm(ln('(/;ooi)) the levels with m = n = 0. 
Employing (11), (12) for m = n = Q one finds 

Vr = jor sin{20), vg — cos{29) (16) 

for the radial and polar components of v, respectively (in the polar coordinates x' — rcos{9), y' = rsm{d)). This 
expression exhibits quadrupolar symmetry. As long as particles are condensing into the state with m = n = I = Q 
they will form the current pattern characterized by (16). This pattern can be thought of as two pairs of vortices of 
opposite vorticity coupled together. 

Above we have shown that no nonadiabatic energy exchange occurs between the RMF and the ideal gas in the trap. 
In the Sec. IV, we will show that the interaction between particles changes this situation. 



III. CONDENSATE CONTAINING LARGE NUMBERS OF PARTICLES IN THE ROTATING FRAME 



In this section we will analyze the case of the condensate containing large numbers of particles in the presence of 
the RMF. The condensate wave function <I> obeys the GGP equation |l9[. For the single-particle Hamiltonian (7) in 
the rotating frame, this equation is 

idt^ = {H^~ ^i)<P + Uom^^ (17) 

where Mq > is the interaction constant and fi stands for the chemical potential. Following the approach 1^, we will 
derive approximate hydrodynamical equations for the condensate in the presence of the RMF. We denote 

$ = x/pe*'^, jd^p = N,, (18) 

where p, <j) and are the density, the phase, and the total number of particles in the condensate, respectively. 
Substituting (18) into (17), one arrives at the expressions 



p - u:{x"d'yp - y'd^^np) + V(pV0) = 0, 

4> ~ oj{x'd'y<t> - y'd.,ncl>) = - Uv^f +^-U-UoP ^ 

2^ 2 



in the rotating frame. The main approximation made in the limit of large Nc is that the term proportional to A(y/p) 
in the second of equations (19) can be neglected ||l^,|8|]. 

In the limit u —>■ oo, one expects to obtain a solution of (19) which is close to that characterizing the TOP JT8|,p[. 
In order to see it, one should separate a rotationally invariant part from the total potential (5). Specifically, 

U = Utop + SU, 

2 2 2 (20) 

Utop = ^i^"' + y") + SU = -'^{x"' - y'^) - x" . 

For the sake of convenience we will omit all primes {x" — > x, y' ^ y) from the coordinates, implying that we are 
working in the frame connected with the RMF unless otherwise stated. Note that Utop is the time averaged potential 
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(TOP) derived in Ref. |£5[. The term 5U describes the deviations of the instantaneous potential (5) from Utop- It is 
not strictly obvious that 5U can be treated as a small correction to Utop- However, the exact results obtained above 
for the single-particle Hamiltonian show that this is true in the limit of large lo at least. Below we will show that if 
Lo » LOo, corrections to the solution (21) caused by 5U remain small for large Nc as well. 
To the zeroth order with respect to 5U , one obtains from Eq.(19) the solution 

0(°) = o, p^°^ = —{^i~UTOp), (21) 

Uo 

which is valid inside the droplet whose radius is determined by the condition p'"' — [p"8P]. We represent p = p^°^ +6p 
where Sp is a small correction due to SU. Correspondingly, we ignore the term W{6pW(l)) in the first equation (19). 
Linearizing Eqs.(19) in Sp, (j) [^i one obtains 

Sp - oj{xdySp - ydjp) + V(p(°) V(/)) = 0, 

(22) 

(f) — Lu{xdy4> — ydx(t>) ~ —UoSp — SU. 

Note that these equations are a close analog to those obtained in ||^ for a trapping oscillator potential which is 
spherically symmetric. The additional feature of Eqs.(22) is the term SU which plays the role of an external force. 
Later we will see that this term under certain conditions can resonantly excite the condensate normal modes with the 
angular momenta i = 1, 2. 

A particular solution corresponding to the symmetry of the driving term SU (20) can be taken in the form 

Sp ^ p'^ix^ - t/2) + 2p'ixv + p{x + p'ly, 

(23) 

= 0^(a;2 - y2) + + c^'^x + ^'(y, 

where pj, p", cj)'/ , I = 1,2 arc the time dependent amplitudes of the dipole {I = 1), and the quadrupolc (I = 2) 
harmonics. Substitution of (21) and (23) into (22) yields 

oj'^ /-I 
(j)l + iluj<j)l + UoPl ^ i-;^) 

,2 2 (24) 
pi + tluipi - I — 0i = 0. 

Un 



for the complex amplitudes 



pi = ip[ + zp[')e^'-\ 0, = (0; + ^0;v-^ (25) 



Note that the dipole amplitudes I — 1 describe essentially the center of mass motion of the whole atomic cloud in the 
trap §. 

If the RMF frequency cu does not change in time, one obtains the steady solutions (p = = 0) 

,^;=0, = p[= f p'l^O, (26) 

and 

This implies that in the limit lo » uJo the corrections due to the RMF to the zeroth order solution (21) [p|Jl^] are 
small. 

Note that Eqs.(27) and (23) indicate that the phase factor exp{ijoXy) , discussed in the Sec. II for the ideal gas situa- 
tion, is not affected much by the interaction as long as ui » uio- Indeed, comparing Eqs.(ll)-(13) with Eqs.(23),(27), 
one finds that the parameter 70 in Eqs.(12), (13), (16) must be replaced by 202 =70 + o{{u;o/^^)^)- When u; cJq, 
the solutions (26), (27) based on the condition p^°^ >> Sp become no longer valid. 

We now consider the case when lu depends on time. For concreteness, we assume that the frequency lu of the RMF 
is modulated as 
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Lu = uJ+ Xsm{Lj't) (28) 

where |a7| >> |A|, and cu' are constants. Accordingly, one finds that the angle between the RMF and the a;-axis in the 
laboratory frame is 

9{t) = ZUt + 4 (1 - cos(w't)) (29) 



Employing the ansatz 



one obtains from Eqs.(24) 



/>,=e-'^(*)^z, p,=e-'^Wp,, (30) 



21 

^i = ^k (31) 

These equations indicate that the resonance condition on w' in (28) is different for the dipole {I = \) and quadrupole 
[l = 2) harmonics. Indeed, given (29) and expanding the r.h.s. of the second equation of (31) in the small quantity 
\\/uj'\ « 1, we get 



e 



lW{t) ^ ilZJt 



ilX 



e'''^«(l + ^(l_coscj'i)). (32) 



UJ 



Then, one obtains that the resonance with the dipole harmonic occurs when the modulating frequency lu' obeys the 
condition 

Uj' = Lo'i =lJ zt LUo- (33) 

Eqs.(31) yield the resonance condition for the quadrupolar harmonic 

= = 2ZU± \/2wo, (34) 

where \/2<^o stands for the frequency of the lowest quadrupolar harmonic of the trapped condensate with large Nc 
Note that ^2 — ~ ti? >> lOq- 
In the following, we will show that the preceding analysis based on the GGP equation does not take into account 
quantum processes of the creation of pairs out of the condensate by the RMF. These lead to forced evaporation 
of the condensate even for zero temperature and steady rotations of the magnetic field. As a consequence, the 
hydrodynamical equations (22) will acquire a dissipative term. 



IV. CENTRIFUGAL INSTABILTY IN THE MANY BODY APPROACH 



In our previous analysis, we neglected quantum fluctuations of the condensate. These fluctuations in the conven- 
tional stable condensate can be thought of as the virtual creation and absorption of pairs. In this regard we note that 
the spectrum (10) has no lower bound, so that the condensate could be unstable with respect to the real creation of 
pairs even though the RMF is steady. Correspondingly, the GGP equation (17) can acquire a dissipative part. 

Gonsider first the case lu ^ oo. The term proportional to lo in Eq.(14) is a direct consequence of the Galiliean 
transformation into the rotating frame. Indeed, the limit uj ^ oo in Eqs.(10)-(13) insures that the term {n — m) is 
the projection of the angular momentum L on the z-axis, so that the w-dependent part in (14) is exactly the Coriolis 
contribution —loL. This implies that no instability should develop because the absence of the lower bound for the 
spectrum is purely a frame of reference effect. Nevertheless the condensate can be considered as being potentially 
able to gain high values of L. In this regard we can employ the rotating frame reasoning ^ (see also 0, Gh.6) 
for the vortex creation in the rotating vessel containing a superfluid. In the frame connected to the vessel rotating 
with the frequency u) around its axis, the vortex energy is — loL^ where Ey and stand for the vortex energy 
in the laboratory frame and the vortex angular momentum, respectively. The vortex can be created spontaneously if 
the Goriolis energy exceeds Ey. However, this argument does not indicate what is the probability for developing this 
centrifugal instability. In fact, in the case of the perfectly symmetric vessel this probability is essentially zero. To 
make the vortex creation real, the vessel must have some irregularities on the walls breaking the rotational symmetry 
so that the angular momentum of the vessel could be transferred to the vortex (or vortices). 
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Returning to our case, we can see that in the case w ^ oo the eigenfunctions (11), (12) of the trap |l^] are 
approaching those of the effective time averaged Hamiltonian (the TOP ||l^) which is axially symmetric. Therefore, 
no centrifugal instabihty of the condensate is expected to occur in this hmit. In other words, no energy exchange 
between the RMF and the condensate happens in the hmit = oo. 

For finite uj, the functions (11), (12) are not eigenfunctions of the operator L. Consequently, the difference n — m 
can no longer be interpreted as the eigenvalue of L. Accordingly the effective vessel can be thought of as having a 
symmetry breaking deformation, which in turn implies that the energy and the angular momentum can now be given 
up to the pairs leaving the condensate into the highly excited states whose energies are tko » fiojo (in physical units). 
In this regard one should distinguish two cases: 1) Tilo » fi and 2) hu! < fi. In the case 1) the pair escapes into 
states lying far from those effectively involved in the formation of the interacting condensate. Accordingly, the pair 
escape process can be treated as an incoherent step in the condensate evaporation. In contrast, in case 2) the escape 
states with the energies « Tiuj are to be renormalized strongly because of the presence of the condensate. This implies 
that the multi-pair processes become significant. Correspondingly, the centrifugal instability should be interpreted as 
a coherent process of vortex formation. In this paper we will not analyze this case. 

The process of the escape of pairs represents the nonresonant quantum evaporation of the condensate induced by 
the RMF. We emphasize the crucial role of the interatomic interaction for realization of this centrifugal evaporation. 
From the point of view of the rotating observer this process can be described as follows: two atoms in the condensate 
interact with each other. As a consequence, they jump to a new pair of single-particle states characterized by large 
quantum numbers, so that their total energy is conserved. Correspondingly, the rotating observer interprets this event 
as a nearly elastic escape of the pair from the condensate. Note that if the eigenfunctions (11) were eigenfunctions of 
angular momentum, there would be a selection rule requiring that the angular momentum of the interacting pair not 
change in the transition. Correspondingly, referring to Eq.(14), one sees that no instability would occur. In fact this 
is not the case for finite to and instability could occur. 

To describe the centrifugal instability effect, we proceed to derive a damping term in the GGP equation for the 
condensate wave function The many body Hamiltonian in the rotating frame is 

J d:K[¥{H^ - /i)* + !^*t^t^^]^ (35) 

where primes are omitted from the coordinates and the Bose operators ^' obey the usual Bose commutation rule. 
The Heisenberg equation is 

idt'^^{H^-fi)<f + Uo^''<f^. (36) 
Taking into account the explicit form (7) for H^, one finds from (36) the current conservation condition 

9t(*t^r) + VJ = 0, (37) 
where the current operator J in the rotating frame is defined as 

J= J-[vl-t(V-A)vI'], 

(38) 

Ax = -iujy, Ay = iLux, Ay^ = 0. 

In the presence of the condensate, the condensate wave function <f> =< 4* >. The noncondensate part = ^I* — $. 
From Eq.(36), one finds [|o|-|| 

5t$ = (iJ^-M)* + "o|$|'$+ 

(39) 

-t-Mo[$* < > +2<I> < VE-'T^-' > + < ^-'T^r'^' >], 

and 

idt^' = {R^ - Ai)*' + «„[$*(«''*'- < >) + 

(40) 

+2<I>(*'%'— < vj/'tvj/' >) + 2$*<i>vE'' + <i)2vl/'t -I- xj/'tvf/'vj/'- < vf'tvi/'vi/' >] 

The condensate wave function $ is normally viewed as an external classical field in Eq. (40) for the noncondensate 
part. Employing the Keldysh technique ||2^, the system (39), (40) can be expressed in terms of the joint dynamics 
of $ and the population numbers of the excitations. In general, this procedure is very complicated |0] (see also the 
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generalized density functional approach However, under certain conditions it becomes possible to eliminate 

the averages from Eq.(39). Specifically, we will make several assumptions and approximations: a) the population 
numbers of the excited states are zero; b) the pairs escaping from the condensate due to the centrifugal effect escape 
from the trap as well; c) the terms leading to powers higher than third in $ and $* in the effective GGP equation are 
omitted; d) in Eq.(40) only the terms linear in ^' and ^''^ are retained in accordance with the Bogolubov approximation 

MM- 

The assumption a) excludes the normal component from the analysis. The assumption b) insures that no normal 
component is building up in the highly excited levels due to the centrifugal escape of the pairs. Note that b) is 
reasonable for high u; and in the presence of the radio-frequency scalpel which provides the evaporative cooling (see, 
e.g., in ||25|,^). Given a) and b), we avoid the necessity to analyze the dynamics of the normal component. Finally, 
from Eqs.(39), (40) under a)-d) we obtain the GGP equation with the dissipation term included (see Appendix B) 

.a* = (H„ - + - /.,'/<ix'G"»(x.,x'i')*nxV,, (41, 

where G*^'"'' (xi, x't') is defined in (B7). Note that if the total number of atoms in the trap were conserved, it would 
not be possible to obtain the dissipation term in (41) in closed form ||]. Below it will be shown explicitly that the 
last term in Eq.(41) would have been zero if either the single particle Hamiltonian 77^ conserved angular momentum 
or if the single particle excitation spectrum were positively defined. 

Multiplying (41) by $ and adding the complex conjugate of the resulting expression, one obtains the generalized 
current conservation condition 



p + divJ ~ —S, 



S 



,2 



dt' [ dx'G('-'2(xt, x'tXxt)p(x't')exp[-2z(0(xt) - 0(x't'))] + c.c, ^^^^ 



where the representation (18) is employed. This equation corresponds to the first equation of Eqs.(22) modified to 
include the dissipation caused by the centrifugal evaporation. 
Integration of Eq.(42) over the whole space yields 

Nc = -^ul J rft'dxdx'Re[G('^)2(xt,x'i')]p(xt)p(x't')exp[-2i(0(xi) - 0(x'<'))] (43) 

where we have used the second relation in Eq.(18). Note that in this equation the integral depends on both the 
density p as well as the phase cj) of the condensate. Accordingly, the coherence of the condensate could be tested by 
analyzing the rate of the quantum vaporization induced by the RMF. Elsewhere, we will consider this possibility in 
greater detail. Presently, let us calculate the quasistatic decay rate assuming that Nc is not large so that one can 
employ the ideal gas ansatz $ = ^/Ncipooo- Note however that the applicability of Eqs.(41)- (43) is not limited by the 
requirement of small {a^JuJ^Nc < 1 |l^, where a is defined in (B5)). When the condition (B5) is satisfied, case 1) 
holds, and one can employ the variational approach jisj for calculating $. 
Making use of Eqs.(ll), (12), we find 



Nc{t) = ~2ul J dt'K{t - t')Nc{t)Nc{t') (44) 

where 

Kit-t')^9{t^t')Y,\Mi2\^Ree-''^'^+'^^^'-''\ M,, ^ /dxV^i^a^o'o (45) 

1,2 

and the summation is performed over the final states of the escaping pair. 

As will be seen below, in the limit uj —^ oo the escape rate is much smaller than the typical time scale in the trap 
corresponding to Uo- Consequently, one can employ the quasistatic approximation that the time dependence of Nc is 
slow. We set Nc{t') — Nc{t) in (44) and after the time- integration rewrite it as 



where 



Nc = -xNl (46) 



X = 27TulY.\Mi2\^6{Ei2). (47) 

[1],[2] 
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In Eq.(47), the notation E12 = Si +£2 for the escaping pair energy is introduced and 6{^) stands for the 5- function. 
This expression accounts quantitatively for the centrifugal effect discussed above. One can see that the condition 
Ml 2 7^ and E12 — 0, where explicitly 

E12 = LOoimi +1712 + 711+ 712) + ^^^02(^1 + h) - L^ini + ^2 - mi - 777,2) = 0, (48) 

can be satisfied simultaneously because the eigenfunctions (11), (12) utilized in (45) and (47) are not the eigenfunctions 
of the angular momentum operator. In what follows, we will show that in the limit of large lo, the dominant 
contribution in (47) comes from the states with the quantum numbers 

mi ~ 7772 ^ li X I2 ^ ni fa n2 ~ — >> 1 (49) 



o 



corresponding to a pair leaving the condensate into the states characterized by large quantum displacements. Re- 
turning to the laboratory frame, this simply means that two atoms absorb the energy 2huj (in physical units) from 
the RMF so that this energy is approximately equally distributed between them. As a result, the pair is transferred 
to highly excited states. The radio- frequency scalpel 25| is assumed to eventually remove this pair insuring the 



condition a) of zero population of the excited states. 

In order to calculate M12 explicitly, we employ the representation (compare with |2l|], Ch.l5) 

^™,„,,,(xi)^™,„,i,(x2) = 5]5r^"^5r"'5^^''V'pfe,(R)^P'fcW'(r) (50) 

pkq 

where the notations p' = m — p, k' — n ^ k, q' = I — q and 



ab _ v^/ -,^p+p' alb\pl{a + b - p)\ 



2—p%a~p')\{p-p')\[b-p + p')\ ^5^^ 

^ = = m^mY+ 7772 , 77 = 771 + 772, I = h + h 

are employed. The summations in Eqs.(50), (51) run over all integer nonnegative numbers obeying the condition that 
all the numbers under the signs of the factorial are nonnegative as well. The relations (50), (51) were derived from 
the explicit representation (11), (12) for the eigenfunctions. Employing Eqs.(50), (51) in Eqs.(45),(47) we obtain for 
Eq.(47) 

2 

X = ^ 51 IV'mn/(0)n7/'000(0)P(5(c^o(777 + 77) + UJ^zl - ^{n ~ 777)). (52) 
rani 

Note that x is exactly zero for the case a = 1 (or a; oo) in Eqs.(ll)-(13). As mentioned above, no escape of the 
pairs occurs in the case when the RMF is so rapidly rotating that the effective trapping (TOP) potential becomes 
axially symmetric. In the limit of large but still finite lo one finds from Eq.(13) 



which implies that the first term only should be kept in (53), and that the exponent e" in Eq.(ll) can be expanded 
in terms of the smallness of lOq/uj. This expansion represents the eigenfunction (11) in terms of the harmonics of 
the angular momentum operator. Keeping the first term only, one finds that each eigenfunction can be effectively 
characterized by three terms: i) the harmonic of the angular momentum operator with the angular momentum {n—m); 
ii) two functions with the momenta (77 ~ 777) ± 2 whose weight is proportional to (ujo/^)^ << 1- This implies that 
in Eq.(47) it is enough to consider the contribution due to the lowest term. Physically this term corresponds to an 
absorption of the energy 2?7Ci; and the angular momentum 2 by a pair of atoms escaping from the condensate. Finally 
one finds 

1 2 2 ^^o^6 (2^)1(777 + 1)(777 + 2) ^^^ 

X = -a uj^ujoz y — ) 2^ Will)^ d{2uJom + 2ujozl -2uj)- (54) 



ml 



A simple analysis shows that most of the contribution to the sum (54) comes from the region of high m, I (see (49)). 
Accordingly, we replace the summation in (54) by integration. Finally, in the chosen limit and chosen units (6) we 
find 
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This expression indicates that the centrifugal escape rate is extremely sensitive to the RMF frequency uj. For the 
parameters employed experimentally in Ref. the estimate of (55) gives a very small number (the corresponding 
lifetime is about lO^s for Nc = 1000), implying that the centrifugal vaporization can be effectively ignored as a cause 
for the condensate escape from the trap. However, with decrease of the ratio w/woj the vaporization rate increases 
strongly. In the case w — > v^Wo from above, the approximation (53) we employed is no longer valid. Accordingly, 
the exact expression (11)-(13) for the eigenfunctions should be utilized in (52). This means that the escaping pairs 
acquire higher (even) angular momenta. As a result, the lifetime of the condensate can become very short. 

Eqs.(39) and (40) can be analyzed for the case of a nonsteady RMF. Especially interesting appears to be the case 
when the RMF excites resonantly the quadrupolar harmonic of the condensate ( see the condition (34)). Generally 
it is natural to expect that this resonance would result in the increase of the vaporization rate as a function of the 
modulating frequency oj' . In the future, we will consider this case in greater detail. 



V. CONCLUSION 



The atomic trap [p^5p] bears features absent in the static traps |^,H,|T6[. These features can be accounted for in 
the frame rotating together with the RMF. For large frequencies of rotation of the RMF, the exact eigenenergies and 
eigenfunctions of the trap [p"5P] approach those characterizing the time averaged potential TOP |l^ having axial 
symmetry with respect to the axis of rotation. For frequencies close to the threshold below which the trapping is 
impossible, the eigenstates loose their axial symmetry and become elongated in the direction perpendicular to the 
RMF (in its frame of reference) . Very close to the threshold a gas of trapped atoms acquires properties of an essentially 
ID system. 

Due to the asymmetry introduced by the RMF, the atom-atom interaction results in the induced evaporation of 
the Bose-Einstein condensate. The time scale for this evaporation is very sensitive to the RMF frequency of rotation. 
For high frequencies, the lifetime of the condensate increases as a large power of lu. Close to the trapping threshold 
the lifetime shortens considerably, implying that the ID gas formed in the trap 1^,^ in this situation is a strongly 
interacting system. 

The RMF can be utilized as a driving force selectively exciting the condensate normal modes. In the limit of large 
numbers of atoms, when Stringari's hydrodynamical approximation is valid, two modes can be excited by the RMF 
whose frequency of rotation is appropriately modulated. The first is a dipole mode accounting for the center of mass 
motion. The second mode which can be excited by the RMF is the lowest quadrupolar harmonic. The resonance 
conditions for the RMF modulation period depend on the averaged RMF frequency, in addition to the eigenfrequencies 
of the harmonics. The effect of the quantum evaporation induced by the RMF opens up a channel for dissipation of 
the condensate normal modes. 
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APPENDIX A: SOLVING THE EIGENPROBLEM IN THE ROTATING FRAME 



The Hamiltonian (9) for the function ip' (x\y' ^ z) — ex.\)(iLijx"y')ip can be written as 

1 

Hi = H^y + = --j^d^ + -^z , 

(Al) 

H.y = -\{dl + dl) + "^^^^y' + \u.'x^ + 2^u.x^y 

where the unimportant constant \/{2uP') and the primes from the coordinates are omitted. The eigenfunctions ^pi{z) 
of Hz are the well known oscillator states. These are represented in (11), (12) by the generating function of the 
Hermite polynomials (see the auxiliary variable ^3 in (11), (12)) so that the ip'{x,y^z) = Lpi{z)il}'j^^{x,y). Performing 
the Fourier transform 
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i,{x,p)= j dye-'PH'{x,y), y = ^ - y , (A2) 

one finds for (Al) 

H.y = -\dl + '-dl + ^^l^p^ + - 2co^/\^J^\xp (A3) 

wliere s — sign(ijj^ — oj'^y). 

In tlie case \uj\ < Woy the Hamiltonian (A3) can be diagonalized by implementation of a real rotation in tlie {x,p)- 
plane. However, no discrete states exist in this case because the effective potential of (A3) turns out to have a saddle 
like shape, with the kinetic part being positively defined. In the opposite limit ( s = 1) the discrete states do exist. 
The diagonalization can be achieved by means of the Lorentz transformation 

X = cosh(t?)^ + sinh(z9)p' 

(A4) 

p = sinh(z9)^ + cosh(i9)p' 

leaving the kinetic part — ^(3^ — dp) invariant. In terms of the new variables {^,p'), (A3) acquires the form 

H.y = + - + ^p'^] (A5) 

where uj± are given in Eq.(lO), and the angle 'd satisfies the equation 

tanh(2^) = ^^^J^^, (A6) 

with rj defined in Eq.(13). The resulting spectrum of the total Hamiltonian is given by Eq.(lO). The eigenfunctions 
of (A5), expressed in terms of the variables, can be converted into the x,p coordinates by means of the relations 
(A4). Finally, performing the inverse of the Fourier as well as the scaling transforms (A2), one finds the normalized 
eigenfunctions (11) - (13). 



APPENDIX B: DERIVATION OF THE GGP EQUATION WITH THE DISSIPATION DUE TO THE RMF 

Under a)-d), Eqs.(39) and (40) simplify considerably. We need to find the lowest order term which contributes to 
the imaginary part of Eq.(39). In Eq.(39), the second term in the brackets does not contribute to the imaginary part, 
so we omit it. The last term in the brackets of Eq.(39) produces the imaginary part. However it can be shown that 
it is proportional to the population numbers of the excited states. Consequently, we omit this term also and rewrite 
(39) as 

idt'^ = {H^ - Aj)$ + Uo|*P* + Uo$*G(xi, xt) (Bl) 



where the equal-time anomalous Green's function ||2^,gl| is defined as Ghit^'t) 'I''(xi)'I''(x'<) > ( the overall 
factor i is omitted here and below in the Green's functions definitions [^0|,^). Note that because the interaction 
potential in (35) is chosen in the (5(x) form, Eq.(Bl) contains G with x = x'. The equation for G can be obtained 
from Eq.(40). In order to accomplish this, we will employ the Bogolubov approximation d). Accordingly, multiplying 
(40) by and taking the average, one finds 

i9tG(xt,x'i) = {H^ - p)^G{^t,x't) + {H^ - /i)x'G(xt,x't) + 

+Mo[2($*(xi)$(xt) + $*(x'i)$(x'i))G(x<,x'i)+ (B2) 

+$2(x'i)G-+(xi, x't) + $2(xi)G+- (xi, x't)], 
where the normal Green's functions for coinciding times are 

G-+(x<,x'i) 1''(xt)*'t(x'i) >, 

(B3) 

G+-{xt,x't) =< *'t(xt)1''(x't) > . 
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In Eq.(B2) the notation (...)x nieans that the single particle Hamiltonian acts on the coordinate x. These 
equations should be supplemented by ones for the normal Green's functions. However, as long as one is only interested 
in deriving the imaginary contribution to (Bl) to lowest order with respect to Uo, significant simplification can be 
achieved. Moreover, for u) » cuo, only the high energy part of the spectrum of the normal excitations contributes to 
the imaginary part of G in Eq.(Bl). Correspondingly, one can neglect the effect of the condensate on this part of the 
spectrum (see the case 1) discussed above). The condition when this assumption is valid can be formulated in terms 
of the smallness of the first term in the square brackets of Eq.(B2) if compared with ujG. In other words, 

n\uj\ » Uo\<^f (B4) 

in physical units. Employing the variational approach jl^ , one can estimate |$p and obtain from (B4) 

aV^A^c < 0.034 f -^j (B5) 

where we have employed the representation Uo — 47ra for the interaction constant Ug in terms of the scattering 
length a in the units (6). Actually, for the parameters of the trap [Q the estimate for Nc gives Nc < 10^. If this 
condition holds, the anomalous Green's function G can be found by iteration with respect to Uo, with the zeroth order 
approximation being zero. Correspondingly, the normal functions (B3) and G should be taken in the zeroth 

order as G^ = and G = (5(x — x') (see conditions a) and b)). Assuming that (B4) (or (B5)) is valid, we finally 
find 

G(xt, xi) = ^iuo J dt' j dx'G^'')^ (xi, x't')$^ (x'i') (B6) 

where the retarded Green's function |^ 

G('-'(xi,xY) = ^(i-i')5]e-'^i(*-*'Vi(x)V^i*(x') (B7) 

1 

is expressed explicitly in terms of the single-particle eigenfunctions (11), (12) and the eigenvalues (10), with the 
summation performed over all the single-particle quantum numbers indicated as 1. In (B7), 0{t) denotes the step 
function. Finally, substitution of Eqs.(B6) and (B7) into Eq.(Bl) yields Eq.(41). 
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